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Abstract 

We investigate the parametric evolution of riddled basins related to synchronization of chaos 
in two coupled piecewise-hnear Lorenz maps. Riddling means that the basin of the synchronized 
attractor is shown to be riddled with holes belonging to another basin in an arbitrarily fine scale, 
what has serious consequences on the predictability of the final state for such a coupled system. We 
found that there are wide parameter intervals for which two piecewise-Hnear Lorenz maps exhibit 
riddled basins (globally or locally), what indicates that there are riddled basins in coupled Lorenz 
equations, as previously suggested by numerical experiments. The use of piecewise-linear maps 
makes it possible to prove rigorously the mathematical requirements for the existence of riddled 
basins. 
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I. INTRODUCTION 



In 1963 Lorenz published a landmark paper which marked the beginning of the modern 
age of nonlinear dynamics [1]. One of the most striking insights of Lorenz was his perception 
that the complicated (in fact, chaotic) dynamics exhibited by the three-dimensional flow for 
some parameter values could be understood in terms of a simpler, low-dimensional discrete 
map. This was obviously a consequence of the strong dissipative character of the equations, 
and it is a widely used procedure in modern nonlinear dynamics research. Lorenz took the 
local maxima of one of the evolving variables of his system and plotted a first return map, 
the so-called Lorenz map. It turned out that the form of this map resembled the so-called 
tent function 1 — 2|a:: — 1/2|) whose properties were already known by mathematicians 

for a long time [2]. 

This resemblance between the Lorenz and tent maps was serendipitous since the latter 
was known to possess ergodic orbits (now called transitive or strongly chaotic [2]), such that 
it provided a sound foundation for the claim that the Lorenz original system would present 
chaotic trajectories in the phase space. Moreover, the correspondence between the butterfly 
attractor and the properties of the Lorenz map furnished an early evidence that chaotic 
dynamics in complex, higher-dimensional systems could be understood by studying simpler, 
low-dimensional systems. 

The Lorenz map is non-smooth thanks to its having a cusp, but it is smooth (albeit 
nonlinear) otherwise. On the other hand, we may have some advantages of using a piecewise- 
linear Lorenz map, also known as the bungalow-tent map [3, 4]. Piecewise- linear maps have 
been intensively studied as models of a variety of physical systems ranging from electrical 
circuits [5] to vibroimpact mechanical oscillators [6]. Moreover, chaotic piecewise-linear 
maps have been shown to be robust in the sense that the chaotic orbits would persist for 
parameter variations, a property not shared by smooth chaotic systems in general [5]. 

In this paper we will explore some dynamical properties of two coupled piecewise-linear 
Lorenz maps. This is a simple system which nonetheless contains some of the fundamen- 
tal properties of coupled map lattices and even oscillator chains. Being piecewise-linear 
such coupled maps have the additional advantage of having a constant Jacobian and thus 
some analytical exact results for, e.g. Lyapunov spectra [7]. This allows one to investigate 
phenomena like complete synchronization of chaos and riddled basins [8, 9]. 
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In brief, a riddled basin is a basin of attraction of an attractor characterized by the 
presence of holes belonging to the basin of another attractor [10, 11]. As a consequence, 
given an initial condition belonging to the basin of some attractor, there will be always some 
points in its neighborhood, no matter how small it is, belonging to the basin of another 
attractor. Hence the physical consequences of riddling are quite serious in terms of our 
ability of predicting what attractor the trajectory originating from a given initial conditions 
asymptotes to. 

Kim and coworkers [12] have suggested the presence of riddled basins in two diffusively 
coupled Lorenz equations, what comprises a six-dimensional system where the synchronized 
orbit lies on a three-dimensional subspace. However, a direct characterization of riddling is 
quite difficult in such high-dimensional systems, calling for indirect ways to evidence riddling. 
One example is to consider the coupling of Lorenz maps, which can give us insight on the 
mechanisms whereby two coupled Lorenz systems can generate riddled basins. However, 
since the original Lorenz map would be hkewise difficult to analyze, a better approach could 
be to use a piecewise linear approximation to it, as we do in this paper. 

Accordingly, our aim in this work is to describe the parametric evolution of riddled basins 
of chaotic synchronization in two coupled piecewise-linear Lorenz maps. We shall describe 
with particular emphasis two situations the parameter ranges for which the coupled maps 
present riddled basins. The parameters to be considered are the slope of the isolated maps 
and the coupling strength, considering two different coupling prescriptions. We show that 
the properties of riddled basins are qualitatively distinct according to the prescription used. 

The rest of the paper is organized as follows: in the second Section we outhne the proper- 
ties of the piecewise-linear Lorenz map. Section HI deals with the synchronization prescrip- 
tions for two coupled maps, with some exact analytical results. The study of riddled basins 
of synchronization are treated in Section IV. The last Section contains our Conclusions. 
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FIG. 1: (a) Lorenz map. (b) Bungalow-tent map for a = 0.45. 



II. THE PIECEWISE-LINEAR LORENZ MAP 



It is a well-known fact that the Lorenz equations, 

dx 

— = a{y-x), (1) 
'^^^ x{(i-z)-y, (2) 

-^ = ^y-i- (3) 

produce phase-space trajectories which converge to a chaotic attractor with an involved 
fractal structure when a = 10^ /3 = 28, and 7 = 8/3 [1]. One way to convey this idea 
is to perform a dimensional reduction by sampling only the local maxima of one of the 
phase-space variables, so as to get a discrete-time variable = maxi=i^{2;(t)}. Then we 
build a first return plot by graphing hn^i versus hn, the result being the so-called Lorenz 
map hn^i = L{hn)j which is a unimodal map with a cusp [Fig. 1(a)]. Such maps display 
a similar dynamical behavior as the tent map {xn^i = 1 — 2\xn — 1/2|) and, in fact, this 
similarity was used by Lorenz to explain the erratic behavior at the chaotic attractor in the 
three-dimensional phase space [1]. 

Another piecewise-linear approximation to the Lorenz map is provided by the so-called 
bungalow-tent map, which has out of four linear segments instead of two, as in the tent 
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function [3, 4]: 



•^n-\-l — fa (-^n) — ^ 



if Xr, G [0,a) 
if Xr, G [a, I) 



a '^^ 

_2a_ I l-3a 
l-2a'^'^ ^ l-2a 

T^(l-:^™) + Hf if^„e [i,l-a) 
i-"^-" ifx„G[l-a,l] 



(4) 



where a G (O, |) is a control parameter. For a = | we obtain the tent map. A linear 
coordinate transformation h = 22x + 32 and the choice a = 0.45 would approximate the 
Lorenz map better than the tent function [Fig. 1(b)]. In the following we will refer to the 
former as the piecewise-linear Lorenz map. 

The two fixed points of the bungalow map are Pi : = and P2 : X2 = I — a. For all 
values of the control parameter a, the fixed points are unstable because their eigenvalues 
satisfy \6i^2 = — ^ 1- Although the map is non-smooth at the fixed point X2j we 

can define the map derivative at the right-hand site of x^^ because this point belongs to the 
interval [1 — a, 1]. 

A useful feature of the bungalow map is that the Lyapunov exponent is a smooth function 
of the control parameter a. In order to obtain this dependence we solve the Frobenius-Perron 
equation 

Pa{x)= / pa{y)6{x- fa{y))dy. (5) 
^0 

to obtain the invariant density [3] 

2a 



Pa (^) = ^ \ X[o,i-a) (x) + -^W-'^^l^-'^^^] ' 



(6) 



where x is the indicator function, i.e., xa {x) = 1 if x ^ and xa (x) = if x ^ A. Thus 
the invariant density is constant in the interval [0, 1 — a) and jumps, at x = 1 — a, to another 
constant value. The Lyapunov exponent is then calculated using 

dfa{x 



A (a) 



/ Pa {X) 

Jo 



In 



1 - a 
2 -3a 



In 



dx 
1-a 
a 



dx. 
1 - 2a 



(7) 



+ 



2 -3a 



In 



2a 



1 - 2a 



In particular, for a = 1/3 we obtain A (1/3) = In 2, which is the expected result for the 
tent map [2], and also turns to be the maximum possible value for this exponent. For other 
values of a, the Lyapunov exponent is always positive, displaying a smooth dependence 
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FIG. 2: Lyapunov exponent of a piecewise-linear Lorenz map versus its parameter a. 

with a [Fig. 2], and assuming the value A = 0.5078 for the approximation to the Lorenz 
map which results from a = 0.45. Hence we shall always be treating a chaotic (actually a 
transitive) map, with no periodic windows irrespective of the value of the parameter a. 

III. SYNCHRONIZATION OF TWO COUPLED PIECEWISE-LINEAR LORENZ 
MAPS 

Coupled maps are paradigmatic examples of higher-dimensional dynamical systems, and 
have been extensively used since the 80s to investigate spatially extended systems where 
there is some kind of interplay between temporal and spatial degrees of freedom. One 
of such phenomena is synchronization of chaos [13]. Coupled chaotic maps can synchronize 
their amplitudes at each time depending on the coupling strength. The synchronized state is 
a low-dimensional attractor embedded in the system phase space, and its stabihty properties 
have been studied with respect to perturbations along transversal directions. Most of these 
investigations, however, are chiefly numerical since the system Jacobian depends on the phase 
space variables. Hence much effort has been devoted to coupled piecewise linear maps, for 
which the Jacobian matrix has constant entries, what simplifies enormously the calculations 
and shed light on the mechanism whereby the synchronized state loses transversal stability 
[9]. 

The simplest problem belonging to this category consists of two coupled piecewise-linear 

6 




Lorenz maps. Remembering that the Lorenz system describes a temporal dynamics for a 
fixed spatial location (the spatial structure being assumed a priori as given), the system of 
coupled maps could be regarded as a rough model of what happens if we consider the time 
evolution at two different spatial locations in the forced convection experiment from which 
we derive the Lorenz equations. In this case a natural choice for the coupling function is 
the diffusive one, i.e. it depends on the difference between the variables for both maps. 

As a matter of fact, Kim et al [12] have considered two identical Lorenz systems diffu- 
sively coupled in their ^-component, or 

^ = a{yi-Xi), (8) 

^ = x,{(3-z,)-y,^ (9) 
dz' 

= XiVi - -iZi + e{zi^x - Zi) (10) 

at 

with z = 1,2, such that z^ = Zi and so on, e being the coupling strength. This system was 
found to present both a synchronized and an anti-synchronized state, defined, respectively, 

by 

{xi = X2, yi = 2/2, zi = Z2), {xi = -X2, vi = -2/2, zi = 22), 

which are contained in three-dimensional invariant subspaces of the six-dimensional phase 
space of the coupled system. There have been found numerical evidences of riddled basins 
for values of e large enough, but a systematic investigation is quite difficult in view of the 
high dimensionality of the system. Hence we can resort to the one-dimensional reduction 
performed by taking the Lorenz map and its piecewise-linear counterpart. 

There are two different forms of coupling two piecewise-linear Lorenz maps, what renders 
two distinct two-dimensional maps, to be referred to as F and G given, respectively, by 

Version 1: ( "'"^^ 1 = F ( ""^ ) = f ["""l + ^ ~ "'^^'l (11) 



Xn+1 \ 


= f{ 


fx. 


Kyn+lJ 










( X, 


.Vn+J 







Version 2: ( ^"^^ ^ ^i-n , f fa + S {y^ x ))\ ^^^^ 

\yn+lj \ynj \Ja[yn + £ [Xn - yn))J 

where 6 and e are coupling strengths which can be different when the coupling is asymmetric, 
and even vanish for unidirectional coupling (a master-slave configuration). As it will turn 
out, the dynamical behavior of the system depends essentially on their sum d = 6 ^ e. 

The basic difference between these kinds of maps is that F does not render the square 
[0, 1]^ invariant under its action, unless both 6 and e are zero, whereas G keeps the square 



invariant when both < 5 < 1 and < £ < 1. In fact, in order to always have F and 
G as well-defined transformations of M^, we must extend the domain of the map variables 
to the whole of M by omitting the restrictions < x and x < 1 in (4). In any case, the 
segment S = {(^, y) | < x = y < 1} is kept invariant and is a subspace where the two 
coupled maps x^^i = fa (x^) and i/n^i = fa {Vn) are synchronized. It will be referred to 
as the synchronization subspace from now on. The orbit in the synchronization subspace is 
obviously the same as a chaotic orbit from the uncoupled map, hence the fixed points in S 
are also Pi = and P2 = 1 — a, located at the intervals 1: [0, a) and 4: [1 — a, 1], respectively. 

In order to check this stabihty of the synchronization subspace we have to investigate 
whether or not a trajectory that starts slightly off S converges to a trajectory in S. This can 
be done by considering the maximum Lyapunov exponent along the direction transversal to 
the synchronization subspace, denoted as A_l. If it is negative (positive) the synchronized 
state is stable (unstable) with respect to infinitesimal transversal perturbations. 

This is only a necessary condition, though, since there may be transversely unstable 
periodic orbits embedded in the chaotic orbit which lies on the synchronization subspace. 
Such transversely unstable orbits may exist even when S is transversely stable (i.e. with 
A_L < 0), and their overall effect is to drive trajectories away from 5, spoiling the convergence 
to the synchronized state. Hence a sufficient condition for a stable synchronized orbit is 
that all unstable periodic orbits in the synchronized state be transversely stable. This last 
requirement can be verified by computing the eigenvalue of the Jacobian matrix along the 
transversal directions such that all of them have moduli less than the unity. 

A. Version I 

Let us now apply these general arguments to the particular case of the coupled piecewise- 
hnear Lorenz maps. Introducing the variables u = x — y and v = x^y there results that the 
synchronization manifold S is given by = 0, and the transversal direction is parameterized 
by the 7;-coordinate. We first consider the coupled system in its version I [Eq. (11)], for 
which the Jacobian matrix is 




(13) 
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where 

'if^ i{x&[0,a), 
^ ^ I if :c G [a, I) , 

' T^a if^e [i,l-a), 
^ ^ ifa^e [l-a,l]. 

The eigenvalues of the Jacobian matrix are given by 

/i2 = c- d, 

with corresponding eigenvectors 




The first eigenvector hes in the synchronization subspace 5, whereas the second eigenvec- 
tor is the transversal one. It should be remarked that, thanks to its being piecewise-linear, 
the coupled map (11) has out of four different results for the eigenvalues, each of them at 
a different interval, according to the value of c given by Eq. (14). Since the elements of 
the Jacobian matrix are constant, their eigenvalues/eigenvectors do not depend explicitely 
on the value of x or y and thus take on the same value for the fixed points Pi = and 
P2 = 1 — a in the synchronized state, belonging to intervals 1 and 4 of [0, 1] and identified 
by c = (1 — a)/a and c = — (1 — a)/a, respectively. 

In Figure 3(a) we plot the logarithm of the transversal eigenvalue (|/^2|) as a function of 
the coupling parameter d for the fixed points Pi and P2 of the coupled maps (11). Hence such 
periodic orbits are transversely stable (unstable) if these numbers are negative (positive). 
We see that, regardless of the value taken on by d, there will be always some fixed point 
with transversal eigenvalue having modulus greater than the unity, although the eigenvalue 
modulus of a given fixed point could be less than unity for some values of d. For the fixed 
point Pi = this interval is di < d < d2^ where 

and for P2 ~ 1 — a, — < d < —di. For a = 0.45, as in Fig. 3 , these intervals are 
h = (-2.23,-0.22) and I2 = (0.22,2.23) for the fixed points Pi and P2, respectively. We 
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FIG. 3: (a) Logarithm of the transversal eigenvalue and Lyapunov exponent of the coupled 
piecewise-hnear Lorenz maps (version I) as a function of the coupling strength d = 5 + e for 
the fixed points Pi = (/i2^^) and P2 = 1 — a {[^^2^)- The intervals for riddling are denoted as 
thicker lines, (b) Magnification of a portion of the diagram. 

conclude that, no matter what the value of the coupling strength d, there will be always 
some transversely unstable fixed point embedded in the (chaotic) synchronized state. 

Note that this analysis has been performed for fixed points only, what it is sufficient for 
the characterization of riddling in this coupling prescription. Other periodic orbits can be 
transversely unstable, however, even when some fixed point is transversely stable. More- 
over, for the intervals Ii and /2, we observe from Fig. 3(a) that one of the fixed points 
is transversely stable, whereas the other one is transversely unstable. This phenomenon, 
called unstable dimension variability, has important consequences in terms of the global 
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dynamics of the system [14, 15], such as the validity of computer-generated trajectories of 
some strongly non-hyperbolic chaotic systems [16] . Moreover, it has applications in chaotic 
systems with more than one attractor, like the ecological problem of two-species competition 
with extinction of either one of the species [17] 



B. Version II 



Now we turn to the coupled maps in its version II, for which the Jacobian matrix is 

(19) 




where c is given by (14), and whose eigenvectors as well as the eigenvalue jii are the same 
as in version I, the transversal eigenvalue /i2 being given by 

l^2 = c{l-d). (20) 

The logarithm of the transversal eigenvalue modulus is plotted in Figure 4(a) as a function 
of the coupling parameter d for the fixed points Pi and P2 of the coupled maps (12). There is 
a common interval for which both fixed points are simultaneously transversely stable, given 
hy ds < d < d^^ where 



and which turns to be the interval Is = (0.18, 1.83) when a = 0.45, as shown in Fig. 4(a). In 
this case, unlike the previous coupling prescription, it is necessary to analyze also period-2 
orbits, since we found that these can be transversely unstable even when the fixed points 
themselves are transversely stable. 
The period-2 points is given by 

^ ^ 3 - 4a ^ ' 2^' ^ ' 

P^'^ ■■ 4^ = fa{x?^) e [1 - a, 1], (23) 

with transversal eigenvalue given by 
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FIG. 4: Logarithm of the transversal eigenvalue and Lyapunov exponent of the coupled piecewise- 
linear Lorenz maps with a — 0.45 (version II) as a function of the coupling strength d — 6 -\- e ioY 
(a) fixed points, and (b) period-2 orbit in the synchronization subspace. The intervals for riddling 
are denoted as thicker lines. 

such that the intervals for which these points are transversely unstable are 

d<l- ^>l+t/4^- (25) 
V2(l-a)' V2(l-«) 

The log-eigenvalue of the period-2 orbit is depicted in Fig. 4(b) as a function of d, showing 
that, for the following intervals of d: [0.4,0.7] and [1.2, 1.5], there will be some transversely 
unstable orbit (either a fixed point or a 2-cycle). 
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IV. RIDDLING AND CHAOTIC SYNCHRONIZATION 



In this Section we shall describe one of the consequences of the synchronized state being 
transversely stable as a whole, but possessing transversely unstable periodic orbits embedded 
on it. For this we focus on the basin of synchronization, which is the set of initial conditions 
in the phase-space which produce orbits which asymptote to the synchronized state as the 
time goes to infinity. It may well happen that this basin of synchronization be riddled, i.e. 
densely punctured with holes (in the measure-theoretical sense) belonging to the basin of 
another attractor. 

The dynamical consequences of riddling are quite severe, for every initial condition we 
may think of are surrounded by an uncertainty ball of smah but nonzero radius, representing 
the observational errors or parametric noise intrinsic to the determination of a state in the 
phase space. Hence, if the basin is densely riddled with regions belonging to another basin, 
such regions will intercept the uncertainty ball, turning uncertain the determination of the 
final state of the system. Even worse, this problem is not going to diminish reducing the 
radius of the uncertainty ball, since riddling occurs at every scale, no matter how small. 

The conditions for the occurrence of riddled basins in a dynamical system are [18]: 

1. There exists an invariant subspace whose dimension be less than of the phase space; 

2. The dynamics on the invariant subspace has a chaotic attractor A; 

3. There is another attractor B, chaotic or not, and not belonging to the invariant sub- 
space; 

4. The attractor A is transversely stable in the phase space, i.e. for typical orbits on the 
attractor the Lyapunov exponents for infinitesimal perturbations along the directions 
transversal to the invariant subspace M are all negative; 

5. A set of unstable periodic orbits embedded in the chaotic attractor A is transversely 
unstable. 

The fact that the synchronization subspace S of the coupled maps is invariant under the 
system dynamics already fulfills condition 1 for riddling. Moreover, if the coupled identical 
systems undergo chaotic motion, their synchronized state is also chaotic, i.e. there is a 
chaotic attractor A embedded in the synchronization subspace (condition 2). For condition 
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3 to be satisfied there must exist another attractor B off the synchronized manifold 5 {global 
riddling). There are many situations, however, for which there is no such attractor, but the 
trajectories off-5 spend an arbitrarily large time wandering through the available phase space 
region before returning to S. In this cases {local riddling) we can still think of this non- 
synchronized transient regime as a second possible behavior and assign to it the properties 
formerly possessed by the second attractor B, namely conditions 4 and 5 of riddling. Hence, 
strictly speaking, local riddling means that only the local conditions for riddling have been 
verified for the system in analysis. 

In practice, however, while most of the conditions for riddling rely on a local analysis in 
the vicinity of the synchronized state, the determination of a second attractor is usually an 
information of a global nature. Except for some cases in which this can be rigorously proved, 
the existence of a second attractor is suggested by numerical simulations of the dynamical 
system. Hence, as far as the local analysis is concerned, we are faced more often with local 
riddling than global one in complex systems. 

A formal way to understand local riddling consists on choosing a hnear neighborhood Sa 
of the synchronization subspace 5, where A ^ 0, which is represented by a strip of width 
2A surrounding it in both sides. If we denote with white pixels those initial conditions 
generating orbits that escape rapidly from and by black points those orbits which do not 
escape and are attracted to the synchronization subspace, we can have the same tongue-like 
basin structure of the globally riddled case, but now restricted to the linear neighborhood 

In the case of global riddhng, it should be mentioned that condition 4 for riddhng implies 
that "most" points, in the traditional two-dimensional Lebesgue measure sense, stay within 
this Sa neighborhood and are eventually attracted to A, as proved in Ref. [10]. This is 
thus compatible with the existence of a positive Lebesgue measure set of orbits in the linear 
neighborhood Sa that move away from it. 

A. Version I 

In order to verify condition 4 for riddling, we have to investigate for which values of 
the coupling coefficient c? = e + 5 the transversal Lyapunov exponent is negative. Let us 
first consider version I of the coupled maps. The Lyapunov exponents along the direction 



14 



-6 -4 -2 2 4 6 

d 



FIG. 5: In the parameter plane a versus d, the black regions have negative transversal Lyapunov 
exponent the coupled piecewise- linear Lorenz maps for a — 0.45 (version I). 

parallel to the synchronization subspace is the same as that for an uncoupled piecewise-linear 
Lorenz map, Eq. (7). The exponent along the transversal direction, on the other hand, can 
be obtained by applying Birkhoff ergodic theorem as 

A_L = lim — > In (^n)l = / \ii\ii2{x) \ p ix) dx 



1 
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1 - 2a 
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1 - 2a 
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2a 



2a 



2 - 3a 



In 



1 - 2a 
-(1-a) 



d 



(26) 



where we have used the invariant density given by Eq. (6). 

In Fig. 3(a) we plot the transversal Lyapunov exponent as a function of the coupling 
strength sum = 5 + £ for a = 0.45. There are two intervals of d for which the transversal 
Lyapunov exponent is negative, thus fulfilling condition 4 for riddling. The boundaries 
of these intervals are given by the values of d for which A_l = 0, yielding transcendental 
equations in view of (26). The regions, in the parameter plane a x d, characterized by 
negative transversal exponents, are depicted in Fig. 5. In the piecewise-linear Lorenz case 
(a = 0.45) such intervals are (-1.224,-1.220) and (0.689, 1.707). 

In the previous Section we concluded that, in version I of the coupled maps, there will 
be always transversely unstable fixed points in the synchronization subspace [see Fig. 3]. 
Besides the fixed point itself being transversely unstable, so are all its pre-images, which 
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form a denumerable infinite set of points, thus condition 5 for riddfing is always fulfilled for 
any value of c?, such that the only remaining condition to be verified is the negativeness of 
A_L, the intervals being thus indicated by Fig. 5 as black regions. 

However, we must remark that condition 5 is of a local nature, whereas the existence 
of a second attractor is an information from the global dynamics of the coupled system. 
Unless we can rigorously prove that there is a second attractor at infinity (or elsewhere) 
for a basin with a nonzero Lebesgue measure the set of mathematical conditions can be 
used to define local riddling. For version I, this proof can be given by writing the coupled 
maps in terms of the variables u = x — y and v = x ^ y. Since the map is defined in 
four intervals, so as the coupled maps in version 1. Let us consider the interval < x < a. 
The dynamics along the transversal direction is piecewise linear, with a slope given by (14). 
Hence, if c? < dc2 = (1 — a)/a the unstable manifold of the unstable fixed point extends to 
infinity, and the pre-images of the fixed point as well. A trajectory that closely follows this 
unstable manifold wih asymptote to infinity, which is the second attractor characterizing 
global riddling. Moreover, this slope must be positive, otherwise the trajectories would 
approach the line x = y but asymptote to infinity as well. This leads to another condition, 
namely d > dci = (1 — 2a) /a. 

For global riddling the basin of synchronization is the set of initial conditions which 
asymptote to the synchronized attractor, with a similar definition for the basin of infinity. 
If we have local riddling, however, it is not possible to define such basins, hence it is more 
convenient to consider the synchronization time r, i.e. the time it takes for an initial 
condition (;ro,yo) to reach the synchronized state (i.e. to approach its vicinity with an 
accuracy of 10~^. 

In Figure 6 we plot in a gray-scale the values of the synchronization time r for the initial 
condition {xq^ yo) placed at the center of the corresponding pixel, for a value of d which fulfills 
the local requirements for riddling. The darker pixels indicate initial conditions whose orbits 
achieve the synchronized state in less than ten iterations, whereas for the lighter pixels this 
time is higher. White pixels indicate orbits which go to infinity. 
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FIG. 6: Synchronization time for initial conditions in the phase plane for a = 0.45 and d = 1.0 
(version I). 

B. Version II 

Now for the version II, the transversal Lyapunov exponent is given by 
1 



A, 



2 - 3a 
= K(a) + ln|l -d\. 



(l_a)ln(i-^j + (l-2a)ln(^" 



2a 



+ ln 1 - d\ 



(27) 



where is a constant for a given a. Figure 4 displays the dependence of A_l with d for 
a = 0.45, showing an interval (o?i,<i2) = (0.442, 1.558) for which the transversal Lyapunov 
exponent is negative, thus verifying condition 4 for riddling. In fact, the boundaries of this 
interval are the roots of A_l = 0, namely 



di,2 = 1 T e" 



-K{a) 



(28) 



straddling the point c? = 1 at which the Lyapunov exponent is minus infinity. 

We have seen in the previous Section that, for the intervals [0.4,0.7] and [1.2, 1.5], we 
can warrant the existence of transversely unstable periodic orbits embedded in the chaotic 
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synchronized attractor, either fixed points or period-2 points, so as to fulfill condition 5 
for riddling. The union of these intervals, on its hand, contain the interval {di^d2) which 
verifies condition 4 for riddhng, and hence the latter ultimately determines the existence of 
riddling. The regions with A_l < are indicated by black pixels in Fig. 7 as a function of 
both a and d. If, besides having A_l < 0, the system also presents riddling, we indicate by 
gray pixels this situation. We observe that, for a = 1/3 (i.e. the tent map case), there is no 
riddling interval, since the system is hyperbolic and A_l < always implies synchronization 
[19]. For a < 1/3, riddling can be explained from the transversal instability of the fixed 
points, whereas for a > 1/3 we have to analyze a period-2 orbit to verify condition 5 for 
riddling. 

In version II of coupling, it turns out that the unit square [0, 1]^ is kept invariant under the 
dynamics of the coupled system if both e and 6 belong to the interval [0, 1]. These constraints 
being respected, there is no attractor in the phase place other than the synchronized one, 
and riddling (if any) is certainly of a local nature. However, if either e or 5 take on values 
outside the interval [0, 1], the unit square is no longer invariant and orbits can, in principle, 
asymptote to other attractor. We have been able to find numerically a second attractor at 
infinity, hence we can actually have global riddling. 

Following the same procedure used in version I, we also plot in Fig. 8 the synchronization 
time for initial conditions and a value of d which is inside the riddling interval. The darker 
pixels, representing "fast" initial conditions, belong to the basin of the synchronized attrac- 
tor, whereas the brighter pixels represent "slow" initial conditions that either spend a long 
time until asymptoting to the synchronized attractor or they may escape to the attractor at 
infinity. The example of riddling given in Fig. 8 shows the regions of "slow" and "fast" initial 
conditions densely intermixed in fine scales. Moreover, the bands of "slow" initial conditions 
have the characteristic tongue-like shape that we find in riddling bifurcation scenarios [20]. 

V. CONCLUSIONS 

The Lorenz equations are a paradigmatic example of deterministic chaos, which has 
stimulated research over the last four decades. One noteworthy feature of coupled Lorenz 
equations is the presumed existence of riddled basins for some values of the coupling strength 
between them, in the sense that the basin structure can be so involved that near a point 
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d 



FIG. 7: In the parameter plane a versus d the black regions have negative transversal Lyapunov 
exponent, whereas the gray regions have in addition a transversely unstable orbit, thus charac- 
terizing riddling. We considered two coupled piecewise- linear Lorenz maps for a = 0.45 (version 

II). 

belonging to some basin there exists a point belonging to other basin, no matter how near 
they are. This is clearly a troublesome feature since any initial condition can be assigned 
with a small but nonetheless finite uncertainty, and the latter can be so large as to turn the 
final state of the system completely uncertain if the basins are riddled. This phenomenon is 
difiicult to analyze in the original Lorenz equations due to the high dimensionality, and thus 
we have considered a low-dimensional model which may shed some light on the dynamical 
behavior of the high-dimensional system, particularly on the existence of riddled basins. 

In this spirit, we analyzed two schemes of coupling piecewise-linear versions of the Lorenz 
map, since exact analytical results can be given for them, e.g. the natural measure of 
chaotic orbits and the corresponding Lyapunov exponents. We found that, as a general 
trend, riddled basins are possible for relatively wide intervals of the coupling strength for 
both versions, what gives us more confidence that the suggestive evidences of riddling in the 
original coupled Lorenz systems are in fact true. There are some differences between both 
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FIG. 8: Synchronization time for initial conditions in the phase plane for a — 0.45 and d — 1.1 
(version II). 

schemes, though. 

When we add the coupling term after applying the map function (version I), the dynamics 
do not keep invariant the unit square. In this case we found that either the synchronized 
state has a (globally) riddled structure or it is not transversely stable and it is not an 
attractor itself. If we couple the variables before applying the map function (version II), the 
dynamics maintains the unit square invariant for some values of the coupling constants. In 
those cases we are sure that there cannot exist a second attractor, yielding local riddling. 
For other values, however, we have not been able to detect this second attractor as well. 
Unlike the previous version of coupling, we found that there are some intervals for which 
there are transversely unstable fixed points and/or period-two orbits. Hence we can have 
presence or absence of local riddling, depending on the parameter values. 

We expect that for coupled Lorenz systems, one would observe a scenario more akin to 
version II, since such simulations typically adopt a diffusive type of coupling. Moreover, the 
more complex dynamical behavior of Lorenz equations makes it possible to assign one or 
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more coexisting attractors, apart from the synchronized state. While bounded attractors 
would be less often probable, the existence of an attractor at infinity is quite common, 
such as in coupled Chua circuits. In those cases our analysis could be extended to the 
characterization of global riddling, as we have done in Version I. 
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